Edge and bulk electron states in a quasi-one-dimensional metal in a magnetic field: 

The semi-infinite Wannier-Stark ladder 



Victor M. Yakovenkod and Hsi-Shcng Goan 
Department of P hysics and Center for Superconductivity, University of Maryland, College Park, MD 20742-4111 
(|cond-mat /9706274|; v.l June 26, 1997; v.2 March 9, 1998; v.3 June 12, 1998) 



We study edge and bulk open-orbit electron states in a quasi-one-dimensional (Q1D) metal subject 
to a magnetic field. For both types of the states, the energy spectrum near the Fermi energy consists 
of two terms. One term has a continuous dependence on the momentum along the chains, whereas 
the other term is quantized discretely. The discrete energy spectrum is mathematically equivalent to 
the Wannier-Stark energy ladder of a semi-infinite ID lattice in an effective electric field. We solve 
the latter problem analytically in the semiclassical approximation and by numerical diagonalization. 
We show explicitly that equilibrium electric currents vanish both at the edges and in the bulk, so 
no orbital magnetization is expected in a Q1D metal in a magnetic field. 

PACS numbers: 71.70.Di, 71.70.Ej, 73.20.At, 73.20.Dx, 76.40.+b 



I. INTRODUCTION 

In a strong magnetic field, the quasi-one-dimensional 
(Q1D) organic conductors of the (TMTSF) 2 X family @] 
(also known as the Bechgaard salts) exhibit very interest- 
ing phenomena, such as magnetic oscillations, magnetic- 
field-induced spin-density wave (FISDW), and the quan- 
tum Hall effect (QHE) (see, for example, review Q). Be- 
cause the Fermi surfaces of Q1D metals are open, these 
phenomena have different mechanisms in Q1D conduc- 
tors compared to more conventional materials with closed 
Fermi surfaces. For example, the QHE exists only in the 
FISDW state, but not in the metallic state of Q1D con- 
ductors Thus far, the theory of the Bechgaard salts 
focused mostly on the bulk electron properties (see, for 
example, review Q). Only recently the edge aspects of 
the QHE in the Bechgaard salts attracted attention |||| . 
An explicit picture of the QHE in the FISDW state in 
terms of the edge states was developed in Ref. ||| . How- 
ever, that work did not take into account possible defor- 
mations of the electron wave functions near the edges. 
In the current paper, we present a detailed study of the 
electron wave functions and energies near the edge of a 
Q1D conductor in the metallic (not FISDW) state. This 
work may serve as a starting point for a more accurate 
theory of the edge states in the FISDW state and their 
role in the QHE. Proper description of the edge states is 
also important for the theory of the cyclotron resonance 
in Q1D metals @. 

The edge states of electrons in a Q1D metal in a 
magnetic field were studied semiclassically by Azbel and 
ChaikinjUH and numerically by Osada and Miura [fiof . 
In Ref. S the WKB quantization condition was applied 
to the problem inconsistently, which resulted in a wrong 
conclusion that the edge states have a discrete energy 
spectrum, whereas the bulk states have a continuous one. 
This statement was also repeated in Ref. |Tl]| . It was 
claimed in Refs. [raM] that the electron edge states pro- 



duce thermodynamic oscillations of magnetization in a 
Q1D metal with an open Fermi surface. In the present 
paper, we clear up the confusion and show that the en- 
ergy of either a bulk or an edge state is a sum of two 
terms, one of which has a continuous spectrum and the 
other discrete (in the approximation where the longitu- 
dinal electron dispersion law is linearized near the Fermi 
energy). The WKB quantization condition determines 
the discrete energy terms of both the edge and the bulk 
states. In Appendix, we explicitly point out a mathemat- 
ical error in Ref. M that led to the wrong conclusions. 



In Sec. IV, we show explicitly that the equilibrium elec- 
tric currents vanish both at the edges and in the bulk, 
so no orbital magnetization is expected in a magnetic 
field. This result is in agreement with independence of 
the bulk internal energy of a Q1D metal with an open 
Fermi surface on a magnetic field 

The Schrodinger equation that we solve analytically 
(Sec. [ij) and numerically (Sec. ill) in order to find the 
discrete part of the electron energy is mathematically 
equivalent to the equations that describe the Wannier- 
Stark ladder |l3| of a semi-infinite ID lattice in a uni- 
form electric field. An analytical solution of this prob- 
lem in terms of special functions was obtained in Refs. 
|p^ , ^5| , but our WKB solution is more general. Our 
results might be useful for interpreting experiments on 
finite-size GaAs-Gai^Al^As superlattices in an electric 
field fill. 



II. ANALYTICAL SOLUTION 

We model the Bechgaard salts by a 2D system that 
consists of ID chains parallel to the x axis and spaced 
at a distance b, their coordinates being y = nb, where 
n is an integer number. The Fermi surface of ID elec- 
tron motion along the chains consists of the two Fermi 
points characterized by the Fermi momenta ±Pp. The 
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energy dispersion law of the longitudinal electron mo- 
tion can be linearized in the vicinity of the Fermi energy: 
en = ±vpp x , where vf is the Fermi velocity, the energy 
en is counted from the Fermi energy, and the longitudi- 
nal momenta p x are counted from ±Pp for the two Fermi 
points: p x — P x T Pf- In this paper, we consider only 
the electron states in the vicinity of the +Pf Fermi point. 
The formulas for the —Pf electrons can be obtained by 
changing the sign of vp. The chains are coupled in the 
y direction by the electron tunneling amplitude t. The 
magnetic field H is applied in the z direction. Choosing 



the Landau gauge, A x — —Hy and A y 



0, we in- 



troduce the magnetic field into the Hamiltonian via the 
substitution p x — > p x — eA x /c, where e is the electron 
charge and c is the speed of light. An energy eigenfunc- 
tion of electron has the factorized form: 



^ 1( tf(a,n)=c*- e /' i ^(n). 



(1) 



The eigenfunctions of transverse motion, </>m(«), are la- 
beled by the discrete quantum number M and obey the 
following ID discrete Schrodinger equation: 

nVt4> M {n) - t[(f>M(n + 1) + <Pm(ti - 1)] = E M <ft M (n), 

(2) 

where f2 is the characteristic energy of the magnetic field: 

n = ebHv F /c. (3) 

Eq. (||) also describes a ID lattice in the uniform electric 
field —Hvf/c in the y direction. This electric field would 
appear in the reference frame moving with the Fermi ve- 
locity vf due to the Lorentz transformation of the mag- 
netic field H. The energy e(p x ,M) of eigenfunction (|l|) 
is the sum of the longitudinal and transverse terms: 



e(p x , M) = v F p x + E M - 



(4) 



We assume that H is not too strong: £1 < 2t. The op- 
posite case Q > 2t, easily treated by perturbation theory 
in the small parameter 2t/Q, requires unrealistically high 
magnetic fields in the Bechgaard salts. 

We consider a crystal that is infinite in the x direc- 
tion and semi-infinite in the positive y direction. The 
wave functions </)j\/(n) are defined at n > 1 with the free 
boundary condition at n = 1. As one can see from Eq. 
(|J), this formulation is equivalent to considering <j>M(n) 
at both positive and negative n with the zero boundary 
condition at n — 0: 



<M0) = o. 



(5) 



We closely follow Ref. [|17| in our treatment of the prob- 
lem. To solve Eq. (|^), we express <^jvf(n) in terms of its 
Fourier transform yjvf (fe): 



<t>M(n) 



k <pM(k) 



dk 
2tt' 



(0) 




FIG. 1. (a) The bulk (curve a) and the edge (curve 6) 
classical trajectories of electrons in the phase space (n,fc*). 
The coordinate n is confined between (Em ± t)/Q for the 
bulk state and between (-Em + t)/Q and for the edge state, 
(b) Solid lines: Classical probability distributions for the two 
trajectories shown in panel (a). Dots: Quantum probability 
distributions |(^m(w)| 2 of the two wave functions correspond- 
ing to the two trajectories. 



Eq. (||) defines the function </>m (n) of the continuous vari- 
able n, which has physical meaning only at the integer 
positive points. The integration in Eq. (|6|) proceeds along 
a certain contour in the complex plane of k. Eq. (^|) sat- 
isfies Eq. (||) provided ifM(k) vanishes at the ends of the 
contour and obeys the following equation: 



iCldifM(k)/dk = (Em + 2t cosfc) ifM(k). 
Solution of Eq. (0) is 

<PM(k)=exp -i £(k',E M )dk' 
Jo 

= exp[— i(Euk + 2isinfc)/£l], 
where the function 

t(k,E M ) = [E M - e ± (k)]/n 



(7) 

(8) 
(9) 

(10) 



is defined for a general transverse dispersion law e±(k), 
whereas Eq. (^ is specific to e±(k) = —2t cosfc. 

When t > Q, integral (§) with tp M (k) from Eq. (§) 
can be taken by the method of steepest descent in the 
vicinity of the points fc* where the derivative in k of the 
phase of the integrand vanishes: 



£(k*,E M ) = (E M + 2tcosK)/Q.. 



(11) 



Eq. (|ll|) can be interpreted as the classical conservation 
law of the kinetic, —2t cosfc,, and potential, nfi, energies 
of electron. If the coordinate n belongs to the classi- 
cally allowed region [(Em — 2t)/Cl, (E M + 2t)/Q], then k* 
is real; otherwise, fc» is complex. Real solutions of Eq. 
(pT[) describe classical electron trajectories in the phase 
space (n, fc„). When Em > 2i, the trajectory lies entirely 
within the region n > and does not cross the boundary 
of the crystal at n — (curve a in Fig. |l](a)). When 
— 2t < Em < 2t, the trajectory reaches the edge (curve 
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FIG. 2. Complex plane of k. Thick lines show the contours 
of integration in Eq. (^|) for three different positions of the 
coordinate n. 

b in Fig. 0(a)). These two types of classical trajectories 
correspond to the bulk and the edge quantum states of 
electrons. The classical motion is periodic both for the 
bulk trajectory, because the end points fc* = ±7r cor- 
respond to the same state, and for the edge trajectory, 
because elastic reflection at point Ub reverses the sign 
of fc* and transfers electron back to point kc- Thus, we 
expect the WKB quantization condition to apply in both 
cases: 



£(k,E M )dk = 2tt(M + 7), 



(12) 



where — l<7<0isa constant, and the integral repre- 
sents the phase space area enclosed by the classical tra- 
jectory. For the bulk and the edge trajectories a and b 
in Fig. |l|(a), these areas are shaded vertically and hori- 
zontally. Contrary to Ref. B , we find well-defined WKB 
quantization areas for both the bulk and the edge trajec- 
tories. 

To derive quantization condition (|l2|) for our model 
formally and to find the constant 7, we need to apply the 
boundary conditions properly. Integral (|^) with <^m(&) 
given by Eq. (^) converges only if the ends of the integra- 
tion contour extend to infinity within the shaded areas 
in Fig. ^, where Imsinfc < 0, and <fiM(k) tends to zero 
at infinity. The right boundary condition in real space, 
4>m{ti) — > at n — > +00, is satisfied provided the con- 
tour of integration starts in area / and ends in area II 
in Fig. ^. Indeed, in the classically inaccessible region 
n — ► +00, solutions of Eq. ( |lT| ) are imaginary. One of 
them, kj± = iarccosh(n'fi/2t), where n' = n — Em/Q, is 
represented by point A in Fig. g. The contour of integra- 
tion connects regions / and II by passing through point 
A. Taking integral (^]) in the vicinity of point A along 
the direction of steepest descent, which is parallel to the 
real axis of k in this case, we find: 



(j) M {n) « exp{-n'[ln(n'n/£) - 1]}/V2tto' 



(13) 



which does satisfy the right boundary condition. 

Now let us calculate <f>M (n) in the classically accessible 
region. In this case, solutions of Eq. (O), represented 



by points B and C in Fig. |[ are real: kc = — fcs — 
arccos(n'f2/2i). The contour of integration connects re- 
gions / and II by passing through points B and C: 



4>M(n) 



x^ M (kB)/V2^^/(2t/ny-n'- 



(14) 



The factors exp(=Fz7r/4) appear in Eq. (|l4|), because the 
directions of steepest descent for points B and C are 
at the angles ^fn/A to the real axis of k. The integral 
from ks to kc in Eq. (|l4|) reflects the change of function 
(||) between points B and C. For the edge states, the 
point n = is classically accessible. To satis fy t he left 
boundary condition (|^), the first line in Eq. ( fl4"|) must 
vanish at n = 0. This generates quantization condition 
L2|) with 7 = —1/4 for the edge states: 



s(-£ M /2t) 



arccos(— Em /2t) 



£,{k 7 E M )dk = 2n[M 



(15) 



Substituting Eq. (0) into Eq. @ gives a transcendental 
equation on Em , which has the following explicit solution 
for the states on the very edge with M t/Vt: 



E M = t{-2 + [(37rfi/2t)(Af - l/4)] 2/3 } 



(16) 



The total number of the edge states is n e dgc = 2t/Q. Eqs. 
( |l5| ) and (jll]) are similar to the edge states quantization 
equations for a closed Fermi surface fl§|| . 

In the classically inaccessible region n < (Em — 2t)/£l, 
solutions of Eq. (|ll|), represented by points D and E in 
Fig. |^, are complex: ko,E = T 71 " + i arccosh(— n'Cl/2t). 
The contour of integration connects regions I and II by 
passing through points D and E: 



<t>M(n) = (-i) 1-exp i-i Z(k',E M )dk' 

x exp{— n'[arccosh(— n'ft/2t) + in}} 

x exp[-V«' 2 - (2t/fl) 2 ]/V2^ ^/n' 2 - (2t/fl) 2 . 



(17) 



The integral between kjj and fc# in Eq. ( p7| ) proceeds 
along the horizontal line [—it, tt and the vertical lines 
[&d,— 7r] and [it, Ue] (see Fig. |); however, the inte- 
grals along the vertical lines cancel. To satisfy the left 
boundary condition, Eq. (^) for a semi-infinite crystal 
or <^Af(n) — » at n — > — 00 for an infinite one, the first 
line in Eq. (|17|) must vanish. This generates quantization 
condition (|12[) with 7 = for the bulk states: 



(,(k,E M )dk = 2-kM. 



(18) 



Substituting Eq. ( |l0|) into Eq. (|18|), we recover the 
Wannier-Stark ladder Jl3| for the bulk states energies: 



E M = Mfi. 



(19) 
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When Eq. ( [l9| ) applies, the function tp M (k) in Eq. (|) 
is periodic: </5 M (fc) = ifM(k + 2n), thus integral (||) can 
be taken only from — tt to it, because the integrals along 
the vertical portions of the integration contour cancel. In 
this case, the bulk wave functions are expressed in terms 
of the Bessel functions J of an integer order: ^(n) = 

J n -M(2t/Q). 

In all cases, as follows from Eqs. (6) and @ with the 
contours of integration shown in Fig. 2, the electron wave 
functions 



4>M(n) = / — exp 

Z7T 



dk 



ink - 

EM/n (2t/n) 



E M k + 2t sin/c 
O 



(20) 



are nothing but the Bessel functions of a general order 
n — Em /f2 in the Sommerfeld representation jl9| . The 
quantized value of the energy Em is determined by the 
boundary condition (||): 



J_ Eu/n (2t/Q.) = 0. 



(21) 



The quantization condition in the form (Bl|) was found 
in Ref. [14|. As shown in Appendix, Ref. W would have 
obtained the same Eq. (^), if mathematical errors were 
not made there. 

The electron wave functions can be expressed in terms 
of the Bessel functions ( |20| ) only when e±(k) = — 2icosfc 
in Eqs. ( |l0| ) and (||), which corresponds to the elec- 
tron tunneling between the nearest neighboring chains. 
Proper description of the Bechgaard salts requires to 
take into account higher harmonics of the transverse dis- 
persion law of electron, such as — 2i'cos2fc, which cor- 
responds to the electron tunneling between the next- 
nearest neighboring chains j2(J. The WKB method de- 
scribed in this section is still applicable for an arbitrary 
transverse dispersion law e±Qc), but the wave functions 
are not the Bessel functions any more. 



III. NUMERICAL SOLUTION AND DISCUSSION 

To verify the semiclassical results, we solve Eq. (|J) for a 
finite number of chains 7i max = 150 3> n c dgc = 2t/Cl = 25 
by numerical diagonalization of the Hamiltonian 



H 







( Q -t 
-t 2fl -t 
-t 30 




V 



\ 



-t (n max - -t 
-t n max O J 



(22) 



The quantum probability distributions |</>m(?i)| 2 of two 
eigenfunctions of Hamiltonian ( |22] ) (the dots in Fig. |l|(b)) 
agree with the classical probability distributions (the 
solid lines in Fig. 0(b)) of the corresponding bulk and 
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FIG. 3. (a) Eigenenergies Em of Hamiltonian (|22j) found 
by numerical diagonalization in the case n max = 150 and 
rkdgc = 2t/Q — 25. (b) Electron dispersion law e(p x ,M) 
(kh. Only the branches with M = 1, 6, 11, . . . , 61 are shown. 

edge trajectories shown in Fig. 0(a). The classical prob- 
ability distributions are proportional to the square of Eq. 
( |l4| ) and are equal to 2/Tv y , where v y = 2tsmk lf /Ti oc 
l/y/ (2t) 2 — (n'fl) 2 is the velocity and T = §dn/v y is 
the period of classical motion. The numerically calcu- 
lated eigenvalues Em of Hamiltonian (^) , shown in Fig. 
||(a), agree with the semiclassical energies found from 
Eqs. ©, ©, and © within less than 1%. As Fig. 
||(a) demonstrates, the energy levels are uniformly spaced 
in the bulk (see Eq. (|l9|)) with the energy (||) propor- 
tional to the magnetic field, but the spacing is different 
and not uniform near the edges. In agreement with Eq. 
( |l6| ) , the spacing of the levels near the edges is sublinear 
(Em oc const + M 2 / 3 ), and the extremal energy levels 
with M = 1 and M = n max are displaced relative to the 
linear extrapolation of the bulk law ( JToj ) by the amount 
AE « T2t. 

Transitions between the energy levels Em in an exter- 
nal ac electromagnetic field constitute the cyclotron reso- 
nance. Because the penetration depth in metals is short, 
we expect the energies (|l6|) of the edge states to show 
up in the surface impedance, as in conventional metals 
|18| . The edge states were neglected in the theory of the 
cyclotron resonance in Q1D conductors (t]]. 

The complete, transverse and longitudinal, dispersion 
law (0) is shown in Fig. |3|(b). It consists of discrete 
branches, each having a continuous linear dispersion in 
p x . The Fermi momenta of the branches, 



P { F M) 



-E M /v F 



(23) 



are defined as the points where the energy e(p x ,M) (Q) 
vanishes. The Fermi momenta of the bulk states are 
spaced uniformly with the distance G = Q/vf = ebH/c, 
but the spacing is different and not uniform near the 
edge. This may have important consequences for the 
FISDW state. The FISDW couples the +P F electrons in 
the eigenstate M with the —Pf electrons in the eigen- 
state M-N [§. As long as Eq. (|l|) applies, the FISDW 
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wave vector Q x — 2Pp — NG exactly matches the dif- 
ference between the Fermi momenta of these states and 
opens an energy gap in their spectrum. N branches of the 
+Pf electrons at one edge of the crystal and N branches 
of the —Pf electrons at the other edge remaingapless, 
because they have no partners to couple with Q|. Even 
though these 2N modes are gapless, electric current is 
not dissipated, because the modes are chiral, and the 
Hall conductivity is quantized: a xy = 2Ne 2 /h, where 
h is the Planck constant Q]. However, the wave vector 
Q x — 2Pp — NG does not match the Fermi momenta ( p3| ) 
near the edges, where their spacing is not uniform. Thus, 
gapless electron pockets should exist there and cause dis- 
sipation in the QHE regime. The size of the pockets may 
be reduced if Q x adjusts to the spacing of the edge states. 
The energetics involved in the latter effect requires a sep- 
arate study. 



IV. EQUILIBRIUM CURRENTS AND 
MAGNETIZATION 



Since the transverse eigenfunctions 4>Ai(n) are real, 
they carry no electric current across the chains. The 
current carried by eigenstates (H) along the chains is 



ip x ,M( n ) = e ( ±y F 



eA x (n) \ , 



(24) 



where m e is the electron band mass, and the signs ± re- 
fer to the ±Pp electrons. To find the total current / at 
the chain n, we sum Eq. ( |2i| ) over M and integrate over 
p x with the Fermi distribution function (at zero temper- 
ature) : 



Pf+PI 



J(n) 



2e E / S (±UF)l ^ (n)|5 



(25) 



2e 2 A x {n) 



(M) 



cm,. 



M=i _ J 



-Pf-p\ 



where the factor 2 comes from the spin of electrons. It 
is understood that the wave functions (jAj (n) and <ff^ in) 
should be used when the integration over P x is close to 
+Pf and —Pf, and some interpolating functions should 
be used for the intermediate values of P x . The result 
does not depend on the contributions far from the Fermi 
surface. 

Taking into account that 



(27) 



we find that the integral J_p dP x in Eq. (^5| ) vanishes. 
The only nonzero contribution to this term comes from 



the deviations pp^ 1 ^ ( p3| ) from the ID Fermi momenta 
±P F : 



— £ -EMi(n)\ 2 + E M \<f> M (n)\ 2 . 

M=l 



(28) 



Taking into account that the eigenfunctions <?^(rj) form 
a complete basis of Hamiltonian (E2L) and using the re- 



+ _ 

M 



lations E 
can be rewritten as 



E~ and Eq. (E7J), we find that Eq. 



^(n\H\n) = ~nSl, 



(29) 



where {n\H.\n} are the diagonal matrix elements of Hamil- 
tonian (p2|). 

Using Eq. (p7j), the term (Eq) can be written as 



2eA x (n) 



j -Pf m=1 



(30) 



Taking into account the completeness relation 
Y) M l^tf ( n )\ 2 — 1 an d integrating over P x , we transform 
Eq. (pjj) into 



2eA x (n)vF 2enfl 
cttH ttH 



(31) 



The two terms (|2^) and (^l|) cancel each other, so that 
the total electric current on any chain n is zero: 



J(n) = 0. 



(32) 



Because the current vanishes everywhere including the 
edges, there is no orbital magnetization (and no de Haas- 
van Alphen oscillations proposed in Refs. in a Q1D 
metal in a magnetic field. Experimentally, no magneti- 
zation was found in the Bechgaard salts in the metallic 
state [^lj (unlike in the FISDW state, where energy gaps 
exist in the electron spectrum). 



V. MAGNETIZATION IN THE CASE OF 
QUADRATIC DISPERSION LAW 

In the previous section, we found that orbital magne- 
tization of the system vanishes identically. That is a con- 
sequence of the linearized longitudinal energy dispersion 
law of electrons in our model. However, for a nonlinear 
dispersion law, magnetization is not necessarily zero. We 
can crudely estimate the change in the bulk free energy 
per one electron at zero temperature generated by an ap- 
plied magnetic field, AF, in the following way. AF must 
vanish when H — > and when t — > 0. (When t = 0, 
the magnetic field has no orbital effect on ID uncoupled 
chains.) Because AF does not depend on the signs of 
H and t, it should be quadratic in £1 = ebHvp/c and 
t in the lowest order. To achieve the dimensionality of 
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energy, we need to divide the expression by a power of 
the Fermi energy £p = PfVf/2. In this way, we find 



AF - t 2 n 2 /e] 



(33) 



Magnetization is obtained by differentiating Eq. ( |33| ) in 
H. 

It is difficult to calculate AF explicitly in the case of 
a weak magnetic field: 51 <C £ <C £_f- In the semiclas- 
sical (WKB) approximation, the bulk free energy of a 
QID metal does not depend on the magnetic field, even 
if the longitudinal dispersion law is nonlinear, as long as 
the Fermi surface is open, and the electron energy spec- 
trum is continuous, not quantized 0). This result is 
related to the Bohr-van Leeuwen theorem, which states 
that partition function in classical statistical mechanics 
does not depend on magnetic field. Thus, in order to 
obtain a nonzero AF, it is necessary to go beyond the 
WKB approximation, which is difficult. 

On the other hand, we can easily calculate AF in case 
of a strong magnetic field: t <C ^ <C Sp, although this 
case may not correspond to the Bechgaard salts in realis- 
tic magnetic fields. In this case, the transverse tunneling 
amplitude t can be treated as a small perturbation to 
the energy spectrum. The second order correction to the 
total energy of the system per one electron at zero tem- 
perature due to a perturbation V is given by the following 
expression: 



AjF (2) = J_ y 

6 <*(e a <e 



(a\V\P)((3\V\a) 
e a - Bp 



(34) 



where N e is the total number of electrons, and the sum 
is taken over the energy eigenstates below and above the 
Fermi energy, labeled by the indices a and j3, respec- 
tively. Treating the transverse tunneling amplitude t as 
the perturbation V and taking into account that its ma- 
trix elements change the longitudinal momentum p x by 
±G: (a k ' x ,M'\V\/3k x ,M) = -tS M ',M±iS(k' x - k x ^fG), we 
find that the sum in Eq. ( |34| ) is restricted to an interval 
of the width G in the vicinity of the Fermi momentum: 



SF <?) = f° 



1 



p e J-G 2vr?i E\\{jp x ) - E\\(p x + G) 



(35) 



In Eq. (|35), the factor of 4 accounts for the two Fermi 
points and two spin orientations, and p e = Akp /2tt is 
the electron concentration per one chain. Using the 
quadratic longitudinal dispersion law Eii(p x ) = (Pf + 
p x ) 2 /2m e , where m e is the effective electron mass, we 
find: 

o 



AF (2) = -— / — 



U 2 

p e J_ G 2wh G(P F + p : +G/2) 
_f_ ( P F + G/2 \ 

n n \p F -G/2j 



(36) 



Expanding Eq. ( pq ) in the small parameter G/Pf = 
VL/2ef and keeping the first two nonvanishing terms, we 
find: 



AF< 2 ) = - — 
2e F 



t 2 n 2 

96e F 



(37) 



The first term in Eq. ( p7| ) coincides with the second- 
order correction due to the electron tunneling between 
the chains in the absence of magnetic field. Only this, 
magnetic- field- independent term is obtained, if the lon- 
gitudinal dispersion law is linearized in Eq. ( p5|) . The sec- 
ond term in Eq. ( |37j ) appears due to nonlinearity of the 
dispersion law and reproduces the result of dimensional 
analysis (^3|) up to a numerical factor. Its negative sign 
indicate paramagnetism. However, because the carriers 
in (TMTSF)2X are holes with a negative m e , the orbital 
response would be diamagnetic in these materials. 



VI. CONCLUSIONS 

In conclusion, the energy of either a bulk or an edge 
electron state in a QID metal is the sum of two terms 
( |) , one of which has a continuous spectrum and the other 
discrete. The discrete part of the electron energy is deter- 
mined by the semi-infinite Wannier-Stark equation (|^). 
We have solved the semi-infinite Wannier-Stark problem 
semiclassically and numerically. The WKB quantization 
condition ( [l2| ) of the electron phase space area (Fig. [I](a)) 
determines the energies of the edge states with the con- 
stant 7 = -1/4 @ and the bulk states with 7 = (|l|). 
The energies are spaced uniformly in the bulk, but not 
near the edges (see Fig. ||(a) and Eqs. ( |l9| ) and ([l6|)). 
These results may be important for the cyclotron reso- 
nance and the QHE in the Bechgaard salts, as well as the 
finite-size GaAs-Gai-^AlxAs superlattices in an electric 
field. We have demonstrated explicitly that the equilib- 
rium electric currents vanish both at the edges and in 
the bulk, so no orbital magnetization is expected in a 
QID metal in a magnetic field in the approximation of 
linearized longitudinal energy dispersion law of electrons. 
We have also estimated the magnitude of orbital magne- 
tization for the quadratic longitudinal dispersion law. 



ACKNOWLEDGMENTS 

VMY is grateful to P. M. Chaikin, A. H. MacDonald, 
R. E. Prange, and S. Das Sarma for useful discussions 
and to also M. Ya. Azbel for drawing attention to Ref. 
H . This work was partially supported by the NSF under 
Grant DMR-9417451 and by the Packard Foundation. 



In this 
(8) of Ref. 



An 



jpendix, 



APPENDIX 



we use the notation of Ref. 



Eq. 



6 



~irMJ v (M)Ji- v {M) = sm(nv), 



(38) 



can be simplified by using identity (9.1.15) from Ref. |^2) 

2 sin(7w) 



J V+1 (M)J- V (M) + J u (M)J-u-i(M) 



ttM 



(39) 



Substituting Eq. @ into Eq. (||), we find: 
J„(M)[J_ y+ i(M) + J- v -x(M)} = -J U+1 (M)J- V (M). 

(40) 

Using the recurrence relation (9.1.27) of Ref. |22|: 



J„ + i(M) + J v -i{M) = ^Ju(M), 



in Eq. (|40|), we find: 
"2;y 



71/ 



J V {M) - J V+1 {M) 



= 0. 



(41) 



(42) 



Using the recurrence relation ([ll]) in Eq. ([l2|) again, we 
find: 



J- v (Af)J v - 1 (M)=0. 
Eq. ( f43| ) is satisfied if either 

J_„(M) = 0, 



or 



J,_x(M) = 0. 



(43) 
(44) 
(45) 



Eq. ( p[ ) is the same as our energy quantization con- 
dition Eq. (pl|). Eq. ( p5| ) describes unphysical electron 
states located outside of the crystal (m < 0) and should 
be discarded. The two sets of eigenvalues, ( |44| ) and (45), 
are completely decoupled and do not repel when cross. 
Thus there should be no gaps in Fig. 2 of Ref. || , no the 
fractional interference between the two set of the energy 
levels, and no diamagnetic oscillations. Contrary to the 
explicit transformation given above, the two sets of eigen- 
values come out coupled via a constant A in Eq. (10) of 
Ref. [El . We conclude that there must be an error in the 
"rather boring calculations" mentioned between Eqs. (9) 
and (10) and leading from Eq. (8) to Eq. (10) in Ref. g. 
The conclusions following Eq. (10) in Ref. || are invalid 
because of the error in this equation. 
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